FRACTIONAL EMBEDDINGS AND STOCHASTIC TIME 
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Abstract. As a model problem for the study of chaotic Hamiltonian systems, we look for 
the effects of a long-tail distribution of recurrence times on a fixed Hamiltonian dynamics. 
We follow Stanislavsky's approach of Hamiltonian formalism for fractional systems. We 
prove that his formalism can be retrieved from the fractional embedding theory. We deduce 
that the fractional Hamiltonian systems of Stanislavsky stem from a particular least action 
principle, said causal. In this case, the fractional embedding becomes coherent. 



1. Introduction 

Fractional calculus lias been widely developped for a decade and its efficiency has already 
been proved in various topics such as continuum mechanics (see |6j), chemisty (see |3j|), trans- 
port theory (see [H]), fractional diffusion (see [7]), etc. 

In [H], a link is drawn with chaotic Hamiltonian systems. Because of the appearance of 
fractal structures in phase spaces of nonhyperbolic Hamiltonian systems, fractional dynamics 
may arise in such systems. Zaslavsky then explained [ISl, chap. 12-13] that time takes on 
a fractal structure, meaning that it can be considered as a succession of specific temporal 
intervals. However, further investigations have to be carried out to understand and clarify the 
link between this peculiar temporal comportment and the fractional dynamics. 

A contribution is done in [H]. As a model problem for the effects of a given distribution of 
recurrence times on the underlying Hamiltonian dynamics, we use Stanislavsky's approach for 
his definition of an Hamiltonian formalism for fractional systems. Indeed, this author looks 
for the effects induced by the assumption that the time variable is governed by a particular 
stochastic process on a given Hamiltonian dynamics. This kind of process contains notably 
the case of the algebraic decay of recurrence times that occurs in the study of chaotic Hamil- 
tonian systems (see [19J). He proves, under strong assumptions, that the induced dynamics is 
fractional and that the structure of the new system looks like the classical Hamiltonian one. 
This allows him to give a definition of an Hamiltonian formalism for fractional systems. 

However, an important property of Hamiltonian systems is that they can be obtained by a 
variational principle, called the Hamilton least action principle (see [2j). A natural question 
with respect to Stanislavsky's construction is to know if his definition of fractional Hamiltonian 
system can be derived from a variational principle. 

In this paper, by using the fractional embedding theory developped in [8J, we prove that 
Stanislavky's Hamiltonian formalism for fractional systems coincides with the fractional Hamil- 
tonian formalism induced by the fractional embedding. In particular, this means that Stanis- 
lavsky's fractional Hamiltonian systems can be obtained by a variational principle. Moreover, 
this fractional formalism is coherent, meaning that there exists a commutative diagram for 
the obtention of the fractional equations. 

In section[2]we discuss Stanislavsky's formalism. Section[3]is devoted to the development of 
the fraction embedding theory using the Caputo derivatives. We obtain a causal and coherent 
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embedding by restricting the set of variations underlying the fractional calculus of variations. 
We also prove that the fractional embedding of the usual Hamiltonian formalism resulting 
from the Lagrangian one is coherent. In section HI we prove that the fractional Hamilton- 
ian formalism stemming from the causal fractional embedding coincides with Stanislavsky's 
formalism. We finally discuss open problems in section O 

2. Stanislavky's Hamiltonian formalism for fractional systems 

2.1. Definition of the internal time. Let Ti, T2, ... be nonnegative independant and iden- 
tically distributed variables, with distribution p. We set r(0) = and for n > 1, T(n) = 
Y17=i The Ti represent random temporal intervals. Let {Nt}t>o = max{n > [ T{n) < t} 
be the associated counting process. We suppose that there exists < a < 1 such that 

/'W^^' t^oo, a>0, 0<a<l. (1) 

Therefore the variables Tj belong to the strict domain of attraction of an a-stable distribution. 
Theorem 3.2 of [13] implies: 



Theorem 1. There exists a process {S{t)}t>o and a regularly varying function b with index 
a such that 

{b{cy Nct}t>o =^ {Sit)}t>o, as 00 

FD 

where =^ denotes convergence in distribution of all finite- dimensional marginal distributions. 

The process {S{t)}t>o is a hitting-time process (see [13]) and is also called a first-passage 
time. From [5], the distribution of {5(t)}t>o, denoted pt, verifies 

C\pt]{v) =E[e-'"'^'^] = E^{-vt"), 

where C is the Laplace transform and Ea is the one-parameter Mittag-Leffler function. It 
follows that 

00 

e-'"^ptix) dt = w^-ie"^'"" (2) 



The process {S{t)}t>o is increasing and may play the role of a stochastic time, which is called 
internal time in [18]. The distribution pt(r) represents the probability to be at the internal 
time r on the real time t. Using this new time, Stanislavsky studies Hamiltonian systems 
which evolve according to S{t). 

2.2. Fractional Hamiltonian equations. We consider an Hamiltonian system, with Hamil- 
tonian H{x,p), and associated canonical equations 

^xit) = d2H{x{t),p{t)), 

If (3) 
= -d,H{x{t),p{t))- 

If i is replaced by S{t), how is the dynamics modified? To answer this question, Stanislavsky 
introduces new variables Xa and pa defined by 

Xa{t) =lK[x{S{t))] = pt{T)x{T)dT, 

p^(t) =nP{Sm = !^ Pt{r)p{T)dT. ^''> 
Furthermore, he assumes that 

diH{Xa{t),Pa{t)) = PtiT)diH{x{T),p{T))dT, 
d2HMt),pUt)) = J^pt{T)d2H[x{T),p{T))dT, 
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which leads to 



Theorem 2. Let {x,p) be a solution of ([3|). Then condition JSj) is verified if and only if 
{xa,Pa) defined by (jl} verifies 



where aX't ^^/^ Caputo derivative defined by 



d2H{Xait),Pa{t)), 
-diH{Xa{t),Pa{t)), 



aVff{t) 

Proof. As X verifies ([3]), we have 
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r 1 - a) 



{t-T)-^f'{T)dT. 



Pt{T)d2H{x{T),p{T))dT 

The Laplace transform of this expression gives 



Jo 



H(x,p) dr 







Pt{r)—x{T)dT 



^\Pt]{s)-^x{T)dT 



„a-l 



oo J 

-x{T)dT from © 



dr 



Given that jC[xa]is) = s" C[x]{s"), we have 



Pt{T)d2H{x{T),p{T))dT 



By taking the Laplace image of this relation, we obtain 



x(0). 



C[oVfx^] {s) 



o'DtXait) 

In a similar way, we also have 



Pt{T)d2H{x{T),p{T))dT. 



(6) 



/•oo 

^Vfp^{t) = - / Pt{T)dlH[x{T),p{T))dT, 

Jo 

and the equivalence follows. □ 

For a presentation of the fractional calculus and its applications, see pT] and [TS]. Hence, 
we will say that a fractional system of the form 

0V^p{t) = f2{xit),p{t)), 

is Hamiltonian in the sense of Stanislavsky if there exists a function H{x,p) such that 

fi{x,p) = d2H{x,p), 
f2{x,p) = -diH{x,p). 

We show that the fractional derivative oVf appears as a natural consequence of the struc- 
ture of the internal time S{t). The fractional exponent a is exactly determined by the be- 
haviour ^ of long time intervals. We note that if we had a > 1 in ([T|), the a-stable distribution 
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would be the Gaussian one, we would have ptir) = dr{t) and then S{t) = t. In this case, inter- 
nal time and real time would be the same. Consequently, for a > 1, the associated derivative 
is the classical one. 

3. Fractional embedding of Lagrangian and Hamiltonian systems 

One important property of classical Hamiltonian systems is that they are solutions of a 
variational principle, called the Hamilton least action principle (see [2]). A natural question 
is to know if the fractional Hamiltonian systems defined by Stanislavky can be derived from 
a variational principle. 

Fractional Euler-Lagrange and Hamilton equations has been first derived in [16], in order 
to include frictional forces into a variational principle. In [1], a fractional Euler-Lagrange 
equation is obtained using a fractional least action principle. This formalism includes left and 
right fractional derivatives. The related Hamilton equations are derived in [4]. However, their 
equations are different from those obtained by Stanislavsky. 

Using the fractional embedding theory developped in [8], we prove that Stanislavsky Hamil- 
tonian formalism stems from a fractional variational principle, called causal^ and moreover that 
this construction is coherent. 

We sum up here the general ideas of the fractional embedding theory for the Caputo deriv- 
ative. Similarly to the left one, the right Caputo derivative is defined by 

The left fractional integral is defined by 

a^r/W = T^ [\t-Tr-'f{r)dr, 

and the right one by 

3.1. Fractional embedding of differential operators. Let f = (/i,...,/p) and g = 
{gi,...,gp) be two p-uplets of smooth functions M^+^ — > MJ. Let a, 6 € M with a < b. 
We denote 0(f, g) the differential operator defined by 

0{i, g){xm = ■ §■9^){x{t), ^x{t),t), (7) 

where, for any functions / and g, (/ • g){t) = f{t) ■ g{t), where • means a product component 
by component. 

The fractional embedding of 0(f, g), denoted fQ,(C'(f, g)), is defined by 

£^{0{{,g)){x){t) =j2{U ■ {aVtrg,){x{t),..,{aVtfx{t),t) 

We define the ordinary differential equation associated to 0{{, g) by 

O(f,g)(x)=0. (8) 
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The fractional embedding <SQ,(0(f, g)) of dS]) is defined by 

3.2. Lagrangian systems. Now we consider a Lagrangian system, with smooth Lagrangian 
L{x, V, u) and u E [a, b]. The Lagrangian L can naturally lead to a differential operator of the 
form ©: 

0{l,L){xm=L{x{t),j^x{t),t). 
Now we identify L and 0{1,L). The fractional embedding of L, £a{L), is hence given by 

£a{L){x){t) = L{x{t),aVr X(t),t). 

In Lagrangian mechanics, the action and its minima play a central role. For any mapping 
g, the action of g, denoted A{g) is defined by 

Aig){x) = f g{x){t)dt. 

J a 

For example, with the identification L = L), the action of L is given by 

A{L){x) = L{x{t),j^xit),t)dt, 
and concerning the fractional embedding of L, the associated action is 

Ai£aimx) = f L{x{t),aVfx{t),t) dt. 
J a 

The extremum of the action of a Lagrangian L provides the equation of motion associated: 

Theorem 3. The action A{L) is extremal in x if and only if x satisfies the Euler- Lagrange 
equation, given by 

diL{x{t), j^x{t),t) - j^d2L{x{t), j^x{t),t) = 0. (9) 

This equation is denoted EL{L). 

This procedure should not be modified with fractional derivatives. Indeed, the strict defini- 
tion of the Lagrangian L does not involve any temporal derivative. The dynamics is afterwards 
fixed with the choice of the derivative T> and the relation v{t) = Vx(t). The variational prin- 
ciple providing the Euler-Lagrange equation uses a integration by parts, which remains in the 
fractional case: 

' UPf fit)] g{t)dt = f fit) m git)] dt 

J a 

+ gib)aV;^'-''^ fib) - /(a),p-(^-") gia). 

We introduce the space of variations 

V^ = {he CH[a, b]) 1 aV;^'-''^ hia) = hib) = 0}. 
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For /i € Fq, we have 
Ai£aiL)){x + h) = Ai£a{L))ix)+ 

b 

[diL + tV^ d2L]ix{t), aVf xit),t) h{t) dt + oih), 
which imphes that the differential of A{Sa{L)) in x is given, for any /i £ Vq, by 

dA{£aiL)){x,h)=j [diL+tV^d2L]{x,aVfx,t)hit)dt, 
= {[diL+tV^ d2L]ixi-),aT>r xi-),-),h), 
where {f,g) = f{t)g{t) dt is a scalar product defined on C^([o, 6]). 

If -B C Va, we will say that A{£a{L)) is i?-extremal in x if for all h ^ E, dA{£a{L)){x, h) = 0. 
So we obtain a first Euler-Lagrange equation: 

Theorem 4. A{£ol{L)) is Va-extremal in x if and only if x verifies 

diL{x{t),aVt x{t),t)+tV'^ d2L{x{t),aVf x{t),t) = Q. (10) 

Proof. A{£a{L)) is Fa-extremal in x if and only if for all /i G Vq, we have 

{[diL+tVt d2L]{x{-),aVt x(-), h) = 0. 

This is equivalent to [diL + {D"^ d2L\{x{-) , aVf x{ )) G . We conclude by noticing that 
= Va^ = {0}, where Va is the adherence of in C^([a, h]), equal to C^([a, h]) entirely. □ 

Equation (llOp will be called general fractional Euler-Lagrange equation and will be denoted 
ELg{£a{L)). Contrary to ([H), two operators are involved here. We will now discuss the 
problematic presence of jl?^ . 

3.3. Coherence and causality. Because of the simultaneous presence of the two derivatives, 
the position of x at time t depends on its past positions, through ^Pf , but also on its future 
ones, through fl?^ . The principle of causality is here violated, which seems crippling from 
a physical point of view. Moreover, we note that ^ can be written in the form with 
f = (1, 1) and g = {diL, —d2L). The fractional embedding £a[EL{L)) of ([9|) is therefore 

diL{x{t),aVt x{t),t)-aV^ d2L{x{t),aV^ x{t),t) = 0, 

which shows that ELg{£a{L)) ^ £a{EL{L)): fractional embedding and least action principle 
are not commutative. So we obtain two procedures providing different fractional equations, 
which seems also unsatisfactory. We are facing a Cornelian choice: shall we preserve causality 
or the least action principle? A possible way to solve this problem is to restrict the space 
of variations. We note Va = {h € Va \ aX'f h = —tT^^ h} and Ka = aPf + , defined on 
C^([a, b\). For any f,g€Va, {Ka f, g) = (/, Ka g). We show that Ka is essentially self-adjoint 
and we obtain a new Euler-Lagrange equation: 

Theorem 5. A{£a{L)) is Va-extremal in x if and only if there exists a function g such that 
X verifies 

diL{x{t),aVt x{t),t)-aVf d2L{x{t),aVt x{t),t) = Ka g. 
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Proof. AiSaiL)) is l/a-extremal in x if and only if [diL+tV^ d2L\{x{-) , aVf x{-)) G V^. Given 
that = {Ker Ka)^ = TmKa, A{£a{L)) is extremal if and only if there exists g such that 

diL{x{t),aV^ x{t),t)+tV^ d2L{x{t),aVt x{t),t) = K„ ~g. 

We conclude by setting g{t) = g{t) + d2L{x{t),aDf x{t),t). □ 

Restricting of the space of variations breaks the unicity of the solution. However, among 
those solutions, there is a single one which remains causal (without the operator tT^^ ), for 
g = 0: 

diL{x{t),J)tx{t),t)-aVfd2L{x{t),aVfx{t),t)=Q. (11) 

Equation ([TT]) will be called causal fractional Euler- Lagrange equation., and will be denoted 
EL,{£^{L)). 

Now causality is respected and we have ELc{£a{L)) = £a{EL{L)). In this case, the frac- 
tional embedding is told coherent, in the sense that the following diagram commutes: 



L -fa(L) 



(causal) LAP 



causal FLAP 



idiL - p2L) = — - (9iL - aVf d2L) = 

where (F)LAP states for "(fractional) least action principle". As aT^l — ~t^b ~ 
can say that the least action is also causal in the classical case. 

However, in the fractional case, the physical meaning of is not clear, but it might be 
related to a reversible dynamics of the variations. Furthermore, this underlines the significant 
role of variations in the global dynamics. 

3.4. Fractional Hamiltonian systems based on fractional Lagrangian ones. There 
exists a natural derivation of an Hamiltonian system from a Lagrangian system based on the 
Legendre transformation. We consider an autonomous Lagrangian system, with Lagrangian 
L{x,v), and we suppose that 

Vx, V I— > d2L{x,v) is bijective. (12) 

The (static) momentum associated to the (static) variable x is p = d2L{x,v). So there exists 
a mapping / named Legendre transformation such that v = f{x,p). The Hamiltonian H 
associated to L is defined by 

H{x,p) =pfix,p) - L{x,f{x,p)). (13) 

It implies diH{x,p) = —diL{x,f{x,p)) and d2H{x,p) = f{x,p). 
Let introduce the function 

(p - d2L{x,v) 
diH{x,p) + diL{x, f{x,p)) 
d2H{x,p) - f{x,p) 

The link between Lagrangian and Hamiltonian formalisms is done through the equation 

Flh{x,p,v,w) =0. (14) 
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The momentum p induces a function p{t) = d2L{x{t),v{t)), which can be considered as the 
dynamical momentum. 

For the classical dynamics, (fT^ becomes 

Flh [x{t),p{t), j^x{t), j^p{t)^ = 0, 

i.e. 

p{t) = d2L{x{t),j^x{t) 
diH{x{t),p{t)) = -diL{x{t),j^x{t) 

d2H{x{t),p{t)) = j^x{t). 

Moreover, if x{t) is solution of the Euler-Lagrange equation ©, we obtain the canonical 
equations 

jx{t) = d2H{x{t),p{t)), 
j^pit) = -d,Hix{t),pit)). 
For the fractional case, the fractional embedding jS]) of (fT^ is given by 

p{t) = d2L{x{t),aVfx{t)), 

diH{x{t),p{t)) = -diL{x{t),aVfx{t)), 

d2H{x{t),p{t)) = aVfxit). 

Then the following result states: 
Theorem 6. If x{t) is solution of the causal fractional Euler-Lagrange equation (fTTI) . we have 

aVfx{t) = d2H{x{t),p{t)), 
aVfp{t) = -d,H{x{t),p{t)). 

These are the equations describing the dynamics of a fractional Hamiltonian system derived 
from a Lagrangian formalism. But Hamiltonian systems can also be considered directly as it 
will be seen. 

3.5. Embedded Hamiltonian systems. Now we consider an Hamiltonian system as defined 
in section El with an Hamiltonian H{x,p) and with equations (E]) associated. The fractional 
embedding ([8]) of JS]) is 

aV?x{t) = d2H{x{t),p{t)), 

aVtpit) = -diH{x{t),p{t)). ^''') 

Furthermore, by indroducing the function 

Lh{x,p,v,w) =pv - H{x,p), 
we see that classical Hamiltonian systems are critical points of the action of Lh defined by 

A{Lh){x,p) = LH[x{t),p{t),j^x{t),j^p{t)^ dt. 
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In the fractional case, the action becomes 

A{SULH))ix,p) = [\h {xit),p{t),aVfx{t),aV?pit))dt. 
J a 

Using the causal fractional Euler-Lagrange equation for Lj^, we obtain 

Theorem 7 (Hamiltonian coherence). Let H be an Hamiltonian function. The solutions 
{x{t),p{t)) of the fractional system (fTSll coincide with causal critical points of the action 
A{£a{LH))- More precisely, the following diagram commutes: 



Lh 



(causal) LAP 



d2H 
-diH 



aVfp{t) 



causal FLAP 
= d2H 



-diH 



Proof. The causal fractional Euler-Lagrange equation for Lh is 

-drH{x{t),p{t))-aV1p{t) = 0, 

aVfx{t)-d2H{x{t),p{t)) = 0, 

which is exactly (fTSl) . 



□ 



So we have coherence between the directly embedded equations and the equations obtained 
by a variational principle. But we have also coherence between this section and the previous 
one, i.e. between the fractional Hamiltonian systems resulting from Lagrangian ones and the 
embedded Hamiltonian systems. 

In other words, the equivalent approaches for Hamiltonian systems in the classical case 
remain equivalent in the fractional case if we use causal variational principles. 

Now we will discuss the link between this formalism and Stanislavsky's one. 



4. Compatibility between the two formalisms 

Condition ([5]) means that the partial derivatives of H commutes with E[-(S'(t))]. This 
condition could just seem of technical order and could appear as unrelated to the real dynamics. 
However, by using the fractional embedding, we can precise the underlying dynamical link, in 
the case of natural Lagrangian systems. 

We consider a natural Lagrangian system, i.e. with a Lagrangian L of the form L(x, v) = 

-mv^ — U{x), and the Hamiltonian H{x,p) derived as in section E31 So we have H{x,p) = 

+ U(x), with p = d2L(x,v) = mv. We suppose that (x,p) is solution of the classical 

2m 

Hamiltonian equations ([3]). We define the associated variables Xa and Pa 

Theorem 8. If Xa is solution of the causal fractional Euler-Lagrange equation (fTTI) associated 
to L, then condition ([5]) is verified. 
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Proof. We set Pa(t) = d2L{xa{t),QVf Xa{t)), i.e. Pait) = m qV" Xa{t). Then, from theorem 
dl]), {xa,Pa) is solution of 

oVfXait) = d2H{Xa{t),Pa{t)), 
0Vfpa{t) = -diH{Xa{t),Pa{t)). 



Moreover, we have 



oo 



d 



Pa{t)=mo'DtJ pt{T)x{T) dr = m J pt{T)—x{T)dT 

Pt{T)mv{T) dr = I pt{T)p{T)dT = po,{t). (16) 
Jo 

So we can replace pa by pa in ([TB]) . which concludes the proof. 

□ 

5. Conclusion 

If we consider the temporal evolution variable of a Lagrangian system as a succession of 
random intervals, and if their density has a power-law tail, then the dynamics of this system is 
fractional. The associated equations can be determined through a fractional embedding, based 
on a least action principle. In order to obtain causal and coherent equations, it is necessary 
to restrict the space of variations. This condition might be seen as a way to cancel the 
finalist aspect of the least action principle. Even if it is still unclear, this model of time could 
notably be appropriated for the description of some chaotic Hamiltonian dynamics. Some 
numerical experiments show that distributions of Poincare recurrence times may possess a 
power-law tail (see p3| chap. 11], [10], [9]). Consequently, the time may be decomposed into 
a succession of recurrence times. For long time scale dynamics, the number of intervals is great 
and the new characteristic time clock may become S{t). This new time takes into account the 
peculiar structure of the recurrence times: if the power-law exponent a verifies < a < 1, 
the long time scale dynamics becomes fractional with the same exponent a. This idea of 
stacked dynamics based on two time scales could be linked with [H], where close results are 
obtained. However, because of the Kac lemma (p2j), which states that the mean recurrence 
time is finite, condition ^ may be valid only locally, near some island boundaries, called 
sticky zones. Further investigations have to be carried on to clarify this point. 
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